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Fractional approximation of solutions of 
evolution equations 

Anatoly N. Kochubei, Yuri G. Kondratiev 


Abstract 

We show how to approximate a solution of the first order linear 
evolution equation, together with its possible analytic continuation, 
using a solution of the time-fractional equation of order <5 > 1, where 
5^1 + 0 . 


1 Introduction 

In many applications (see, for example, IZIIH) we encounter evolution equa¬ 
tions, which are too singular for application of the standard semigroup ap¬ 
proach. A possible remedy is to consider the equation on a scale of Banach 
spaces and to use an approach known as Ovsyannikov’s method; see [3] for 
its history and complete references. 

Barkova and Zabreiko [1] extended this method to evolution equations 
with the Caputo-Djrbashian fractional derivative of order a > 0. In 
particular, they showed that the Cauchy problem considered in this setting 
becomes less singular with the growth of a. For example, it can happen 
that a solution u of the hrst order equation u[ = Au exists only on a hnite 
time interval, while the existence of solutions of the equation = Au, 
1 < 5 < 2 can be guaranteed for all f > 0. Therefore such solutions are 
natural means of approximating solutions of the initial hrst order equation. 

In this note we prove an even stronger result - the “fractional approxi¬ 
mations” approach, as 5 1 -I- 0, not only the solution u, but its 
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maximal possible analytic continuation. As we will show, this follows from 
Hardy’s theorem [5] about approximate analytic continuation obtained by 
summing certain divergent series. 


2 Cauchy problems 


Consider the Cauchy problem 

— = Au, m(0) = uq, 


in a scale of Banach spaces a; e [0,1), X'^ c A" for u' < u". Here A is a 
linear operator, which is bounded from X'^ to X" for each couple of indices 
with uj' < O’", and 


For each couple (a;', w") and any initial vector uq e such that Auq e A^/, 
there exists a local solution in Xi^n of the form 



see [21 [7]. 


Let 1 < 5 < 2. Consider also the Cauchy problem with the Caputo- 
Djrbashian fractional derivative 


= Aus, M5(0) = Mo, ^^(O) = 0. 
Under the same assumptions, it has the solution (see DP) 



which exists for alH > 0 (the above series equivalent to the iteration process 
from [T] converges on any hnite interval). 

In fact, the solution u{t) can be continued to a holomorphic function on 
the disk {f e C : |t| < while the function 



71 = 0 
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is extended to an entire fnnction. Both are with values in Note that 


j-n 

usit^/^) = V A„(<5) (A-uo) ^ 

„ n\ 


n=0 


where A„(h) 


r(5n+ !)■ 


3 Analytic continuations 

Definition (see 0). Let f{z) be a holomorphic function on a neighborhood 
of the point z = 0 determined there by a convergent power series. The Mittag- 
Leffler star G{f) of the function f is a domain obtained from C as follows: 
draw a ray from the origin to each singular point of the function f, and cut 
the plane along the part of the ray located after the singular point. 


00 

For example, the star of the function ^ z^ 

n=0 


1 

1 - Z 


is C\[l, oo). 


Theorem. Suppose that the solution u{t) is continued along rays to a single¬ 
valued holomorphic function on G{u). Then 

— *u{t), teG{u), 

uniformly on any closed bounded domain inside G{u). 


Proof. By Theorem 135 (Section 8.10) from [5] (evidently valid also for 
vector-functions), it is sufficient to prove that the function 


Mz) 


^T{5n+lf ’ 

n=0 ^ ' 


(5 > 1 , 


is entire and tends to-uniformly, as h ^ 1 -I- 0, in each closed bounded 

1 — z 

domain not intersecting the semi-axis [1, oo). The entireness property follows 
from the Stirling formula. The above convergence is obvious for \z\ < 1. Due 
to the uniqueness theorem for holomorphic functions, it will be sufficient to 
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prove that tends to some holomorphic function uniformly on compact 

subsets of C\[0,oo). 

The function (ps can be represented as a Wright function: 

] ^ [r(i + k)r z^ 

r(i + 6k) k\ 




with parameter A = h — 2, which means that — 1 < A < 1 (for information 
about the Wright functions see [6], 1.11). 

Let us use the integral representation for the Wright functions. In our 
case we have 




i+ioo 

j_ (- r(s)[r(i-s)]^ 

2711 J r(i — hs) 

i-ioo 


{—z) ^ds. 


Since r(s)r(l — s) = we hnd that 


;+200 


Mz) = 




r(i-s) 


2i J sin(7rs)r(l — Ss) 


(—z) ^ds, I arg(—z)! <(1) 


(see [B], (1.11.21), (1.11.22)). Note that for 6 close to 1 this representation 
becomes valid for any 2 : outside the semi-axis (0, 00 ). 

By a corollary of the Stirling formula ([6], (1.5.14)), for s = I -1- iy, 
\y\ —> GO, we have 


r(i - s) 

r(i-fo) 




5-1 

2 


exp 


i)|y| j 


Noticing also that | sin(7r(4 + iy)\ = \ cos{i7ry)\ = | \e^y + we see the 

possibility to pass to the limit in (3.1), as 5 ^ 1 -I- 0. Thus, for z ^ [0, 00 ), 
uniformly on compact sets. 




l+ioo 




sm TTS 


ds. 


The right-hand side is a holomorphic function on C\[0, 00 ). 


□ 
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